home *** CD-ROM | disk | FTP | other *** search
/ Mac-Source 1994 July / Mac-Source_July_1994.iso / C and C++ / Libraries / stringsearch / bmsource / uf.rev.d12.c < prev    next >
Encoding:
C/C++ Source or Header  |  1991-05-06  |  3.3 KB  |  165 lines  |  [TEXT/MPS ]
  1. /*
  2.     search routine generated by gen.
  3.     skip=uf, match=rev (using revr), shift=d12
  4. */
  5. /*
  6.  * The authors of this software are Andrew Hume and Daniel Sunday.
  7.  * 
  8.  * Copyright (c) 1991 by AT&T and Daniel Sunday.
  9.  * 
  10.  * Permission to use, copy, modify, and distribute this software for any
  11.  * purpose without fee is hereby granted, provided that this entire notice
  12.  * is included in all copies of any software which is or includes a copy
  13.  * or modification of this software and in all copies of the supporting
  14.  * documentation for such software.
  15.  * 
  16.  * THIS SOFTWARE IS BEING PROVIDED "AS IS", WITHOUT ANY EXPRESS OR IMPLIED
  17.  * WARRANTY.  IN PARTICULAR, NEITHER THE AUTHORS NOR AT&T MAKE ANY
  18.  * REPRESENTATION OR WARRANTY OF ANY KIND CONCERNING THE MERCHANTABILITY
  19.  * OF THIS SOFTWARE OR ITS FITNESS FOR ANY PARTICULAR PURPOSE.
  20.  */
  21.  
  22. #ifndef    CHARTYPE
  23. #define    CHARTYPE    unsigned char
  24. #endif
  25. #define    MAXPAT    256
  26.  
  27. #include    "stats.h"
  28.  
  29. #ifndef    TABTYPE
  30. #define    TABTYPE    long
  31. #endif
  32. typedef TABTYPE Tab;
  33.  
  34. static struct
  35. {
  36.     int patlen;
  37.     CHARTYPE pat[MAXPAT];
  38.     Tab delta[256];
  39.     Tab delta1[256];
  40.     Tab delta2[257];
  41. } pat;
  42.  
  43. prep(base, m)
  44.     CHARTYPE *base;
  45.     register m;
  46. {
  47.     CHARTYPE *skipc;
  48.     register CHARTYPE *pe, *p;
  49.     register int j;
  50.     register Tab *d;
  51.     register Tab *d2;
  52.     register q1, tp, t, qp, jp, kp;
  53.     Tab f[256], f1[256];
  54.  
  55.     pat.patlen = m;
  56.     if(m > MAXPAT)
  57.         abort();
  58.     memcpy(pat.pat, base, m);
  59.     skipc = 0;
  60.     stats.len = m;
  61.     d = pat.delta;
  62.     for(j = 0; j < 256; j++)
  63.         d[j] = pat.patlen;
  64.     for(p = pat.pat, pe = p+m-1; p < pe; p++)
  65.         d[*p] = pe-p;
  66.     d[*p] = 0;
  67.     skipc = (CHARTYPE *)p;
  68.     d2 = pat.delta1;
  69.     for(j = 0; j < 256; j++)
  70.         d2[j] = m;
  71.     for(j = 0; j < m; j++)
  72.         d2[base[j]] = m-1-j;
  73.     d2 = pat.delta2;
  74.     for(j = 1; j < m; j++)
  75.         d2[j] = 2*m-j;
  76.     for(j = m, t = m+1; j > 0; j--, t--){
  77.         f[j] = t;
  78.         while((t <= m) && (base[t-1] != base[j-1])){
  79.             if((m-j) < d2[t])
  80.                 d2[t] = m-j;
  81.             t = f[t];
  82.         }
  83.     }
  84.     q1 = t;
  85.     t = m+1-q1;
  86.     qp = 1;
  87.     for(jp = 1, kp = 0; kp < t; jp++, kp++){
  88.         f1[jp] = kp;
  89.         while((kp >= 1) && (base[jp-1] != base[kp-1]))
  90.             kp = f1[kp];
  91.     }
  92.     while(q1 < m){
  93.         for(j = qp; j <= q1; j++)
  94.             if(m+q1-j < d2[j])
  95.                 d2[j] = m+q1-j;
  96.         qp = q1+1;
  97.         q1 += t-f1[t];
  98.         t = f1[t];
  99.     }
  100. /*for(j=1; j<=m; j++)printf("[%d]=%d ", j, d2[j]); printf("\n");/**/
  101.     d2[0] = m+1;        /* the case where the match succeeded */
  102. }
  103.  
  104. exec(base, n)
  105.     CHARTYPE *base;
  106. {
  107.     int nmatch = 0;
  108.     register CHARTYPE *e, *s;
  109.     register Tab *d0 = pat.delta;
  110.     register k;
  111.     register CHARTYPE *p, *q;
  112.     register CHARTYPE *prev = pat.pat+pat.patlen-1;
  113.     register Tab *d2 = pat.delta2+1;
  114.     register k1;
  115.  
  116.     s = base+pat.patlen-1;
  117.     e = base+n;
  118.     memset(e, pat.pat[pat.patlen-1], pat.patlen);
  119.     while(s < e){
  120. #ifdef    STATS
  121.         k = d0[*s];
  122.         stats.jump++;
  123.         while(k){
  124.             stats.jump++; stats.step[k]++;
  125.             k = d0[*(s += k)];
  126.             stats.jump++; stats.step[k]++;
  127.             k = d0[*(s += k)];
  128.             stats.jump++; stats.step[k]++;
  129.             k = d0[*(s += k)];
  130.         }
  131. #else
  132.         k = d0[*s];
  133.         while(k){
  134.             k = d0[*(s += k)];
  135.             k = d0[*(s += k)];
  136.             k = d0[*(s += k)];
  137.         }
  138. #endif
  139.         if(s >= e)
  140.             break;
  141. #ifdef    STATS
  142.         stats.slow++;
  143. #endif
  144. #define    RH    s
  145.         for(p = prev, q = RH; p > pat.pat; ){
  146. #ifdef    STATS
  147.             stats.cmp++;
  148. #endif
  149.             if(*--q != *--p)
  150.                 goto mismatch;
  151.         }
  152.         nmatch++;
  153.     mismatch:
  154.         k = d2[p-pat.pat];
  155.         k1 = d0[*q];
  156.         if(k < k1)
  157.             k = k1;
  158. #ifdef    STATS
  159.         stats.step[q+k-s]++; stats.jump++;
  160. #endif
  161.         s = q+k;
  162.     }
  163.     return(nmatch);
  164. }
  165.